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ABSTRACT 

We report on a 40 ks Chandra observation of the TeV emitting high mass X-ray binary HESS J0632+057 
performed in February 201 1 during a high-state of X-ray and TeV activity. We have used the ACIS-S camera 
in Continues Clocking mode to search for a possible X-ray pulsar in this system. Furthermore, we compare the 
emission of the source during this high state, with its X-ray properties during a low state of emission, caught by 
a 47 ks XMM-Newton observation on September 2007. We did not find any periodic or quasi-periodic signal 
in any of the two observations. We derived an average pulsed fraction 3a upper limit for the presence of a 
periodic signal of ^ 35% and 25% during the low and high emission state, respectively (although this limit 
is strongly dependent on the frequency and the energy band). Using the best X-ray spectra derived to date for 
HESS J0632H-057, we found evidence for a significant spectral change between the low and high X-ray emission 
states, with the absorption value and the photon index varying between Nu ~ 2.1 — 4.3 x lO^^cm"^ and 
r ~ 1.18 — 1.61. At variance with what observed in other TeV binaries, it seems that in this source the higher 
the flux the softer the X-ray spectrum. 

Subject headings: X-rays: binaries — stars: individual (HESS J0632-H057) 



1. INTRODUCTION 

X-ray binary systems are one of the few astronomical ob- 
jects that, under some conditions, are expected to appear as 
point-like TeV emitting sources when observed by instru- 
ments having the current sensitivity. That is the case of the 
three well-established members of the class, PSRB 1259-63 
(Aharonian et al. 2005a), LSIh-61°303 (Albert et al. 2006), 
and LS 5039 (Aharonian et al. 2005b, 2006). Among these 
three systems, PSRB 1259-63 is composed by a 48ms pul- 
sar in a 3.4 ys orbit with a Be companion, LS Ih-6I°303 and 
LS 5039 have a ~26 and 4 days orbital periods, and host a Be 
and an O start, respectively. Unfortunately the nature of the 
compact objects in these two binary systems is still unkown. 

Of the many tens of unidentified sources discovered in the 
Galactic Plane H.E.S.S. survey (the survey done with the 
High Energy Stereoscopic System), HESS J0632-I-057 is one 
of only two unidentified very-high-energy gamma-ray sources 
which appear to be point-like within the experimental resolu- 
tion (the other is coincident with the gravitational centre of 
the Milky Way; Aharonian et al. 2007). 

Follow-up observations of HESS J0632H-057 with 
XMM-Newton have revealed an X-ray source (XMMU 
J063259.3-H054801) coincident with the TeV detection as 
well as with the massive star MWC 148 (spectral type BOpe), 
at a distance of ^ 1.5 kpc (Hinton et al. 2009). The chance 
coincidence of a massive star within the 1 arcsec error box of 
the brightest X-ray source in the XMM-Newton observation 
has been quantified by the latter authors to be of the order 
of 10-^ XMMU J063259.3H-054801 has in addition been 
found to have a hard non-thermal X-ray spectrum, and 
significant variability on hour timescales. These features are 
similar to those found in the many X-ray observations of 
LSlH-6I°303 and LS5039 (SidoU et al. 2006, Esposito et 
al. 2007; Kishishita et al. 2009; Rea et al. 2010, 2011 and 
references therein), strengthening the association between the 
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TeV source and XMMU J063259.3H-054801/MCW 148, and 
hence establishing the possibility for HESS J0632H-057 to be 
the fourth TeV binary. 

Observations conducted with the Giant Metrewave Radio 
Telescope (GMRT) and the Very Large Array (VLA) revealed 
the radio counterpart of HESS J0632H-057: a point-like, low- 
flux, variable radio source at the position of MCW 148 was 
detected in both 1280 MHz with GMRT and 5 GHz with 
VLA (non simultaneously), with an average spectral index of 
-0.6 which can be generated by synchrotron-emitting elec- 
trons (Skilton et al. 2009). 

Subsequent observations of HESS J0632H-057 with the 
VERITAS array detected no significant signal from it, exclud- 
ing that the source is a steady gamma-ray emitter (Acciari et 
al. 2009, Maier et al. 2009). Simultaneous to these TeV 
observations, an X-ray campaign conducted with 5vy//f-XRT 
revealed a significant variability in the X-ray emission (Fal- 
cone et al. 2010). From this early X-ray monitoring there was 
no signature of an orbital periodicity for HESS J0632-H057, 
however, assuming that the spectral variability is due to an or- 
bital modulation, the orbital period was estimated to be larger 
than ^ 54 days. A similar conclusion was reached by Arag- 
ona et al. (2010) through optical spectroscopy of the massive 
star Confirming the previous constraints, Bongiorno et al. 
(201 1) recently reported an orbital period of 320±5 days, us- 
ing data from the continuing long-term Swift-XKT monitor- 
ing campaign. By analogy with LSIh-61°303 and LS5039, 
the compact object hosted in HESS J0632H-057 might then be 
a young non-accreting pulsar, or an accreting compact object 
(black hole or neutron star) driving a jet. 

On 201 1 January 23rd, Swift-XRT detected a rise in the X- 
ray flux of HESS J0632H-057 (Falcone et al. 2011). The flux 
increase was a factor of ~ 3, which appeared similar to the 
rises that occurred ^ 320 and ^ 640 days before. We now 
know that this is the orbital period of the binary (Bongiomo 
et al. 2011). Motivated by this increase in X-ray activity, TeV 
observations were conducted by VERITAS on 2011 Febru- 
ary 7-8, which detected the source at higher TeV flux than 
during the previous VERITAS campaigns (Ong et al. 2011). 
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Fig. 1. — XMM-Newton (left panel) and Chandra (right panel) back-ground 
and 700 s, respectively. 

These results were soon confirmed by the MAGIC collabo- 
ration (Mariotti et al. 2011). Contemporaneous with the X- 
ray increase, further radio observations were also conducted, 
which suggested the presence of slightly extended radio emis- 
sion at milli-arcsec scales coincident with the position of the 
Be companion (Moldon et al. 2011). However, simultane- 
ous optical observations conducted from 201 1 January 5 and 
February 24, revealed no significant change in the radial ve- 
locities of the Be companion star (Casares et al. 201 1). 

During this high TeV/X-ray emission period we have con- 
ducted a 40 ks observation with the Chandra X-ray Obser- 
vatory kindly granted to us using the Discretionary Director 
Time (DDT). Here we report on the spectral and timing char- 
acteristics of HESS J0632-I-057 during this high X-ray emis- 
sion state, and we compare them with an archival XMM- 
Newton observation (Hinton et al. 2009) performed during 
a low X-ray emission state. 

2. OBSERVATIONS AND DATA ANALYSIS 

The Advanced CCD Imaging Spectrometer (ACIS) cam- 
era on board of the Chandra observatory (Weisskopf et al. 
2000) observed HESS J0632H-057, on 201 1-02-13 (start time 
21:15:13 (UT); Obs-ID 13237) for an exposure time of 40ks 
in Continuos Clocking (CC) mode (FAINT). We have cho- 
sen this mode since it provides a time resolution of 2.85 ms, 
suitable for searching for fast pulsations. The CC mode also 
provides imaging along a single direction. The data analy- 
sis mimics that performed in Rea et al. (2010, 2011). Data 
were reprocessed using the CIAO software (ver. 4.3) and the 
Chandra calibration files (CALDB ver. 4.4.3). The source 
was positioned in the back-illuminated ACIS-S3 CCD at the 
nominal target position (RA 06:32:59.30; Dec H-05:48:01.00; 
Hinton et al. 2009). Standard processing of the data was per- 
formed by the Chandra X-ray Center to level 1 and level 2. 

We corrected the times for the variable delay due to the 
spacecraft dithering and telescope flexure, starting from level 
1 data and assuming that all photons were originally detected 
at the target position. We excluded hot pixels, bad columns, 
and possible afterglow events. Finally, photon arrival times 
are in TDB and were referred to the barycenter of the Solar 
System using the JPL-DE405 ephemeris. 

Events in the 0.3-lOkeV energy range were extracted from 
a small region of 5 x 5 pixels around the source position for 
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subtracted light-curves in the 0.3-10, 0.3-2 and 2-10 keV band, binned at 1400 s 

timing analysis, so as to reduce the background contamina- 
tion in the timing analysis. The source spectrum was instead 
extracted from a rectangular region of 5x25 pixels around 
the source position, with the background being taken inde- 
pendently from a source-free region in the same chip. 

Response matrix (RMFs) and ancillary response (ARFs) 
files were produced first creating a weighted image, re- 
binning by a factor of 8. We used this re-binned image to 
build the RMF file using the mkacisrmf tool, with an en- 
ergy grid ranging from 0.3 to 10 keV in 5 eV increments. 
Using this RMF and the aspect histogram created with the 
aspect solution for this observation (asphist), we gener- 
ated the appropriate ARE file for the source position. The 
source ACIS-S count rate in the 0.3-10 keV energy band was 
0.265 ± 0.004 counts s^^ (all errors in the text are reported at 
90% confidence level, unless otherwise specified). 

We have also re-analized the XMM-Newton observation 
reported in Hinton et al. (2009), in search for X-ray pulsa- 
tions, and with the aim of comparing the low and high X-ray 
state of HESS J0632-I-057 using the two best available spectra. 
T\\t XMM-Newton Observatory (J ansen et al. 2001) observed 
HESS J0632H-057 on 2007-09-17 (start time 01:22:50 (UT); 
obs-id 0505200101) for an exposure time of 47 ks with the 
EPIC-pn in Prime Full Extended Windowed mode. The addi- 
tion of the two MOS cameras gave consistent results, we then 
decided to use only the pn camera. This observing mode re- 
sults in a timing resolution of 199.2 ms, a 27.2 x 26.2 arcmin^ 
field of view and only 2.3% of out of time events. Data have 
been processed using S AS version 1 1 .0.0 with the most up- 
dated calibration files (CCF) available at the time the reduc- 
tion was performed (April 201 1). Standard data screening cri- 
teria were applied in the extraction of scientific products. 

The last part of the observation was affected by proton flares 
which we have removed in our final products, except for the 
event file used for the pulsation search (which is not affected 
by such kind of flares), to use as much source photons as pos- 
sible. Cleaning the events from the proton flare results in a 
good exposure time for spectral analysis of 27 ks. 

We have extracted the source photons from a circular re- 
gion of 15" (such to avoid any chip gaps), centered on the 
source point spread function (PSF), and the background from 
a larger region far from the source but in the same CCD 
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Fig. 2. — Left column: Timing analysis of the XMM-Newton data taken during a low intensity state. Right column: same as the left column but for the new 
Chandra data taken during a high intensity state of the source. Top panels: Power spectra of the two observations with the relative 3(t coherent signal detection 
limit. Bottom panels: 3(t limits on the pulsed fraction of a detectable signal in different energy bands. 



(we have re-scaled the background spectrum to take into ac- 
count the different extraction region with respect to the source 
one). For the spectral analysis we used only photons with 
PATTERN<4, and FLAG=0. The source EPIC-pn count rate 
in the 0.3-lOkeV energy band was 0.064 ± 0.002 counts s"^ 
We report in Fig. [T] on the light-curves of the two observa- 
tions we study here, which clearly show also the short-term 
X-ray vai-iability of HESS J0632H-057 (see also Hinton et al. 
2009; Acciari et al. 2009; Falcone et al. 2010; Bongiorno et 
al. 2011) 



3. RESULTS 

3.L Timing analysis: search for pulsations 

We searched for periodic and quasi periodic signals in the 
Chandra and XMM-Newton observations performing a series 
of Fast Fourier Transforms (FFTs; van der Klis 1989). Given 
the length of our two observations (^40 ks each), the tim- 
ing resolution of the Chandra (2.85 ms) and XMM-Newton 
(199.2 ms) observations, and the number of counts of our 
observations, we could search for periodic signals in the 
0.005 — 175 Hz and 0.005 — 8 Hz frequency range, respec- 



tiveljj^ Furthermore, for both data-sets we performed the 
search in the 0.3-lOkeV energy band, and also dividing the 
entire data set in two energy bands (0.3-2 and 2-lOkeV). 

For the Chandra data we performed an average over 7 FFTs 
with a bin time of 2.85 ms (see Fig.l2]right column), resulting 
in about 2,097,152 frequency bins for each of the 7 averaged 
power spectrurrrl The resulting power-spectrum had a dis- 
tribution with iTdegree of freedom (dof). 

For the XMM-Newton data we could perform a single EFT 
with a bin time of 60 ms (over-sampling by a factor of 3 the 
timing resolution of the instrument), resulting in 716,732 fre- 
quency bins in the power spectrum. The resulting power- 
spectrum had a distribution with 2 dof (see Fig.plleft col- 
umn). 

Note that given the long orbital period of HESS J0632H-057 
(Bongiorno et al. 2011) with respect to the exposure time 
of the observations we report here, we do not need to de- 
modulate the photon arrival times for the orbital motion as 
we did for the 4 day TeV binary LS 5039 (Rea et al. 201 1). 

In calculating the 3(t detection upper limits reported in 

' Note that in the search we have oversampled the XMM-Newton timing 
resolution by a factor of 3 to reach the 8 Hz upper bound (see Israel & Stella 
1996). 

^ We could not use less FFTs because of computing limitations. 
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Fig . 3 . — Left panel: Spectral energy distribution of HESS J0632+057 during the liigh (red; Chandra) and the low states (black; XMM— Newton) for an absorbed 
power-law model with all parameters free to vary among the two spectra. Right panel: residuals of the fits of an absorbed power-law with (from top to bottom) 
A'^H and r equal for both spectra, Nn equal and F free, A^h free and F equal, and all parameters free (see text for details). 



Fig. [2] we took into account for the number of bins searched, 
for the different dofs of the noise power distribution, and for 
the red-noise (Vaughan et al. 1994; Israel & Stella 1996; Rea 
et al. 2010). We did not find any periodic or quasi-periodic 
signal in any of the two observations. As a final try we at- 
tempted a joint search using both observations in a single 
power spectrum, with again no detection of any signal. How- 
ever, note that the ^3 years time-span between the two obser- 
vations would likely hamper the detection of a periodic signal 
without the knowledge of its first derivative, and of the system 
precise orbital parameters. 

We computed the 3(t upper limits on the amplitude of a si- 
nusoidal signal (which we define as pulsed fraction (PF)), 
according to Vaughan et al. (1994) and Israel & Stella (1996). 
These limits range in the 0.3-10 keV energy band between 
PF <22-34% (0.005 - 175 Hz) and 32-48% (0.005 - 8 Hz), 
in the high and low X-ray state respectively (see Fig.l2]bottom 
panels). We also infer the same limits as a function of the en- 
ergy band, which given the lower number of counts causes the 
energy-dependent PF limits to be slightly larger than those 
derived using the whole energy range in the two datasets. 

3.2. Spectral analysis 

For the spectral analysis we binned both the Chandra and 
XMM-Newton spectra such as to have at least 50 counts per 
spectral bin (see FigjsJ^ We first fitted both spectra together 
with an absorbed power-law (phabs and powerlaw mod- 
els under the XSPEC version 12.5.0 spectral modeling pro- 
gram). All the absorption values we report here are referred 
to abundances from Anders & Grevesse (1989) and photo- 
electric scattering cross-section from Balucinska-Church & 
McCammon (1998). Fitting the two spectra with the same 
model with all parameters except the normalization value 
equal among the two spectra, results inaxt ~ 1-07 (190 dof; 

' Not e tha t the small difference in the derived spectral parameters we re- 
port in i]3.2|with respect to Hinton et al. (2009), is most probably due to a 



different spectral re-binning used in the analysis. Anyway, all the results of 
our analysis of the XMM-Newton observation are consistent within a 99% 
confidence level with those reported in Hinton et al. (2009). 



with Nil = (4.0 ± 0.2) X 10^1 cm-2, T = 1.56 ± 0.03). Al- 
though the xt value is acceptable, the residuals of the XMM- 
Newton observation are clearly not satisfactory (see Fig[3]). 
This is mirror of the fact that the large difference in counts 
among the two spectra make such a joint spectral fitting (and 
consequently the xt value) to be dominated by the Chandra 
data. 

To solve the problem of having such bad residuals, we 
analize possible spectral variabilities between the two spec- 
tra. In particular, we re-fit the spectra keeping the absorp- 
tion value (A^h) equal among the two spectra and the power- 
law photon index (F) and normalization free to vary. Al- 
though resulting in an acceptable xt = 1 07 (189 dof; with 
A^H = (4.0 ± 0.2) X 10^1 cm-2, Fcxo = 1-57 ± 0.03 and 
Fxmm = 1-46 ± 0.05), the residuals of the XMM-Newton ob- 
servation were again not good (see Figj3j. We then did a fur- 
ther trial keeping equal the F while leaving free the A^h and 
normalization. Again we got an acceptable xt — 1-08 (189 

dof; with A^H,cxo = (4.0 ± 0.2) X lO^l Cm-^, AfH.xmm = 

(3.4 ± 0.3) X 10^1 cm-2, and F = 1.56 ± 0.03), but bad 
residuals. 

Leaving all parameters free to vary for both observations 
we get a good xt ^ 0-98 (188 dof) and better, flat-looking 
residuals for both XMM-Newton and Chandra. We then con- 
sider this as the best spectral modeling for both observa- 
tions. In particular we find the following spectral param- 
eters: for Chandra (during the high X-ray emission state), 
A^H.cxo = (4.3 ± 0.2) X 10^1 cm-2, F^xo = 1-61 ± 0.03, and 
an absorbed (unabsorbed) 0.3-lOkeV flux of (3.2 ± 0.2) x 
lO^^^erg s^^cm^^(4.5 x lO^^^erg s^^cm^^). Foi XMM- 
Newton (during a low X-ray emission state): A^H.xmm = 
(2.1 ± 0.4) X 10^1 cm~2, Fxmm = 1.18 ± 0.08, and an 
absorbed (unabsorbed) 0.3-lOkeV flux of (5.1 ± 0.9) x 
lO^^^erg s^^cm^^(5.7 x 10"-^^erg s^^cm^^). 

To investigate further the spectral change between the two 
observations, we plot the statistical contours for the absorp- 
tion value and the power-law photon index (see Fig|4]l, which 
clearly shows that the spectral variability among the two X- 
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Fig. 4. — Contour plots of the absorption value (A^h) and pboton in- 
dex (F) derived from the fit of the XMM-Newton spectrum taken while 
HESS J0632-H057 was in a low X-ray emitting state, and of the Chandra spec- 
tram during a high X-ray state. Note that the A^h refers to abundances from 
Anders & Grevesse (1989) and photoelectric scattering cross-section from 
Balucinska-Church & McCammon (1998). 

ray emission state is significant at > 99% confidence level. 

Assuming a distance of 1.5 kpc (Hinton et al. 2009), the 
system had a luminosity change between the low and the high 
X-ray emission state of about an order of magnitude, with 



the frequency, energy range and emission state; see \ 2>.\ and 
Fig|2]for details). The limits we derived for the X-ray pulsed 
fraction of HESS J0632H-057 are similar to those derived for 
the archetypical TeV binaries: LS Ih-61°303 (Rea et al. 2010) 
and LS 5039 (Rea et al. 201 1). Note that also the only firmly 
established TeV binary containing a pulsar, PSRB1259-63, 
does not show X-ray pulsations (Chernyakova et al. 2009). 
This result shows that in the pulsar scenario, the pulsar emis- 
sion itself cannot be the main responsible for the X-ray emis- 
sion of these TeV binaries, which is instead likely dominated 
by the wind-wind or intra-wind shock (unless also in this case 
the putative pulsar is never pointing to us during its rotation). 

Also for HESS J0632H-057, we find here that there exists 
a significant spectral variability among the low and high X- 
ray emission states (the flux of the new Chandra observa- 
tion is compatible with that observed by Swift in the same 
period; Bongiorno et al. (2011)). In particular, comparing 
the two best available spectra, which were taken at different 
orbital phases, we can see > Scr variability in the source spec- 
tral parameters. Furthermore, while F increases from the low 
to high state by ^ 40% (hence the spectrum softens accord- 
ingly), the A^H increases by more than a factor of 2. Note 
also that for HESS J0632H-057 we find that the higher the flux, 
the steeper the X-ray spectrum, which is at variance with the 
cases of LS5039 and LSIh-61°303, where the opposite be- 
havior is found (Kishishita et al. 2009; Rea et al. 2010). In 
comparison with the other TeV binaries, we are tempt any- 
way to relate this spectral change with the orbital phase of 
HESS J0632H-057. However our current ignorance on the or- 
bital parameters of this system (beside the 320 days period; 
Bongiorno et al. 201 1) hamper for the moment any classifica- 
tion of the low and high states in terms of periastron, apastron, 
or any of the conjunctions. 



■^low 



1.54 X lO^^ergs-iandLhigh ~ 1-21 x lO^^ergs-i. 



4. CONCLUDING REMARKS 

We reported on a Chandra observation during the 2011 
February increase of X-ray/TeV emission of the new TeV bi- 
nary HESS J0632H-057, which allowed us to perform the first 
detailed X-ray timing and spectral analysis of this source dur- 
ing its high state. As a comparison, we also studied the best X- 
ray data available in the archive for HESS J0632H-057, taken 
while the source was in its low X-ray emission state (previ- 
ously published in Hinton et al. 2009). We do not find any pe- 
riodic or quasi-periodic signal from this system in any of those 
emission states, deriving a 3(7 upper limit on the X-ray pulsed 
fraction of HESS J0632H-057 of -30% (highly dependent on 
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